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• This work is aimed at modeling the channeling phenomenon on a downdraft gasifier. 

. These phenomena have been modeled and compared to experimental data. 

• Results confirmed the capability of this approach to model the channeling phenomenon. 

• The channeling phenomenon depends on the loading method of the gasifier. 
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Downdraft stratified gasifiers seem to be the reactors which are most influenced by loading conditions. 
Moreover, the larger the reactor is, the higher the possibility to stumble across a channeling phenome¬ 
non. This high sensitivity is due to the limited thickness and superficial placement of the flaming pyro¬ 
lysis layer coupled with the necessity to keep all the zones parallel for a correct running of this kind of 
gasifier. This study was aimed at modeling and investigating the channeling phenomenon generated 
by loading condition variations on a 250-kW e nominal power gasification power plant. The experimental 
campaign showed great variations in most of the plant outputs. These phenomena were modeled on two 
modified mathematical models obtained from literature. The results of the models confirmed the capa¬ 
bility of this approach to predict the channeling phenomena and its dependency on the loading method. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Despite the several advantages that characterize the fluidized 
bed gasifiers (Knoef, 2005), fixed bed gasifiers are maintaining a 
high interest in the biomass market due to their simplicity, espe¬ 
cially when operated with air as a gasification agent (Knoef, 
2012; Higman and van der Burgt, 2008). However, fixed bed gasifi¬ 
ers are affected by various problems, such as channeling and fuel 
bridging. 

In the bibliography, channeling can be found related to different 
parameters, such as the residence time (Yamazald et al., 2005), 
superficial velocity (Tarhan, 2003) or devolatilization (Simone 
et al., 2013). In this study channeling was established working on 
the loading frequency of the reactor. It was proved that continuous 
loading of a gasifier reactor avoids strong fluctuations in the gasi¬ 
fication performance (Reed and Das, 1988). These fluctuations are 
common in batch-working gasifiers due to the continuous chang- 
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ing of the reaction conditions (Beedie et al., 1996). Unfortunately, 
continuous feeding does not guarantee the complete absence of 
these oscillations: in concrete terms, most of the applications 
approximate the continuous feeding through a series of small 
batch loads. Every small load causes a pressure or temperature 
drop in the reactor, changing its internal conditions, furthermore 
they can force the air to pass through a few paths, causing 
channeling. 

Gasifier design influences its sensitivity to operating conditions 
(Mechanical Wood Products Branch, 1986; Basu, 2010); the more 
distant the combustion zone is from the loading area, the less 
heavily the combustion reactions are influenced by the loading 
operations. Stratified downdraft gasifiers are characterized by the 
proximity of the flaming pyrolysis zone to the top surface; more¬ 
over, the drier the biomass is entering the reactor, the thinner 
the drying zone becomes (Reed and Das, 1988). This presumes a 
higher sensitivity of stratified downdraft gasifiers performance to 
the variations of loading operations, if compared with other 
typologies. 

This study was aimed at modeling a 800 kW tfI stratified reactor 
under two different loading conditions, one of which was 
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characterized by strong non-homogeneity of the reaction condi¬ 
tions inside the gasifier (due to the channeling effect). The system 
analyzed here consisted of a reactor able to produce up to 
250 kW el , sending the filtered syngas to two IVECO Cursor IC 
engines. 

During the experimental campaign, the plant was run for five 
days, forcing the auger to load the reactor for 15 s every time the 
laser-meter, used for bed height control, gave the low-reactor sig¬ 
nal. After three days of rest, the test was repeated reducing the 
charging time to 4 s. The system was set up fixing the power out¬ 
put to 160 kW c( in order to slow down the loading operations from 
the nominal power conditions. Running the plant at partial power 
output allowed precise monitoring of the following parameters: 
biomass consumption, tar, water and charcoal production, superfi¬ 
cial and inner temperatures, impeller spinning velocities and en¬ 
gine power output. Furthermore, the higher heating values of tar 
and charcoal were evaluated by experimental analysis for both 
operating conditions. 

Results showed that the loading frequency was strongly corre¬ 
lated with the channeling, the explanation of the phenomenon is 
reported in Section 2. 



2. Methods 

A sketch of the analyzed power plant is reported in Fig. 1. It con¬ 
sists of a gasifier reactor fed with poplar woodchips; two diesel en¬ 
gines each set with power output to 80 kW e ; a cylindrical dryer, 
which dried the biomass using the exhaust and the cooling air of 
the engines. Before reaching the engines, the syngas was cleaned 
passing through a cyclone, a water scrubber and a battery of 24 
electrostatic filters. The engines needed 5% of diesel oil to assure 
the ignition of the whole fuel mixture; this fuel fraction could be 
switched from diesel to vegetable oil after a start-up period. The 
reactor is shown in Fig. 2. It had a diameter of 0.92 m and it was 
0.8 m high. In this kind of reactor, a coexistence of different reac¬ 
tion layers takes place (Pedrazzi et al., 2012; Blasi, 2000; Blasi 
and Branca, 2013). A stirrer placed between the flaming pyrolysis 


I UMID BIOMASS 
COOLING AIR f 



Fig. 2. Reaction design. 

layer and the reduction layer supports the biomass mixing and 
avoids the formation of a biomass pile. The observed phenomena 
is explained here after looking locally at the processes that oc¬ 
curred on the reactor surface when the fresh biomass was loaded. 
A large amount of biomass can choke the flaming pyrolysis pre¬ 
venting the correct air flux to the zone. The air entering the reactor 
was diverted into a few paths (channels). 

The previous phenomenon, well known as channeling, was re¬ 
ported in literature (Yamazaki et al., 2005; Tarhan, 2003; Simone 
et al., 2013; Dasappa et al., 2011) and it was implemented in a 
model specifically set up for hypothesizing the gasifier composed 
of two zones with different reaction conditions. The experimental 
data about water, tars and charcoal were compared with the model 
outputs in order to find out the reaction conditions of the two 
zones. Possible solutions capable of reducing the tar production 
through working on loading frequency are suggested here. 

2.1. Experimental campaign 

The experimental campaign discussed here is based on the 
experience acquired during a preliminary experimental investiga¬ 
tion (Allesina et al., 2011 ) where the whole gasifier-filters-engines 
system was outlined an it was found to be characterized by a high 
inertia. The gasifier took hours to heat up properly and the filtering 
system discharged tars produced a long time before. For this rea¬ 
son the experimental campaign was designed, dividing measure¬ 
ments into two different categories: those obtained from the 
steady state plant (i.e. biomass consumption, syngas flow, temper¬ 
atures, impeller spinning velocities, power output) and those ob¬ 
tained from the whole experimental campaign (i.e. water, char 
and tar production). The campaign lasted five days for every load¬ 
ing condition. Every day, 4 h after the plant start-up, the steady 
state measurements were taken. 

2.2. Measurements 

The gas flow rate was measured by a Pitot-tube flowmeter in¬ 
serted into the pipe connecting the filtering system to the engines. 
The measuring point was placed about 30 diameters of rectilinear 
pipe away from the closest curve. In this way the average velocity 
in the pipe ensured a fully developed turbulent flow of syngas; this 
guaranteed good reliability of the measurements. 
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For a precise measurement of the biomass consumption, the 
tests were started after 4 h from plant start-up, while the plant 
was fed from the dryer. During the tests, the dryer was stopped 
and the reactor was fed with a known quantity of biomass previ¬ 
ously dried and weighed. The average mass of poplar wood chips 
used for every steady-state test was about 420 kg. In order to 
maintain the moisture level under control, three samples of bio¬ 
mass were tested for every plant run using a PCE - MB 200 mois¬ 
ture meter. All the tested biomasses had a moisture content of 
under 10%. The water produced in the gasification process was col¬ 
lected from the water scrubber settling tank, and then it was dis¬ 
charged to a 1 m 3 storage tank; after every test the water level 
was measured and the storage tank was weighed. The top readings 
from the reactor temperatures were obtained with three N type 
6 mm thick thermocouples. The thermocouples were placed at 
10, 20 and 30 cm from the top reaction surface, the thermocouple 
tips were placed about 3 cm inside the reactor. The top surface 
temperature in the reactor was monitored continuously with a 
dual laser IR thermometer. Tar and char productions were col¬ 
lected at different points of the filtering system. The charcoal 
was discharged from the bottom of the reactor and the cyclone, 
while the tars were collected from the electrostatic filters or sepa¬ 
rated by floating or sinking in the water scrubber settling tank. 
These processes of separation and collection avoided an instanta¬ 
neous monitoring of waste production. For this reason, the ac¬ 
quired data refers to five days of experimental campaign. The 
instantaneous power output as well as the impeller spinning veloc¬ 
ities were obtained from the plant management software. 


2.3. Ashes analysis 


The biomass used in this study consisted of poplar wood chips. 
The main properties of the biomass are reported in Table 1. During 
the maintenance operation that followed the last test, some 
slugged ashes were found in the reactor. For this reason a sample 
of biomass was heat-treated in a lain for 3 h at 550 °C in order to 
reduce it to ashes. Then a parallelepiped of pressed ashes was 
tested in an optical vertical dilatometer. The dilatometer was 
heated up to 500 °C and then the ashes sample was added (flash 
test mode as suggested by Paganelli and Sighinolfi, 2009). The tem¬ 
perature in the instrument rose to 850 °C at a rate of 80 °C/min and 
then the heaters were turned off, letting the sample cool down. 
Fig. 5 reports three different trends: the temperature of the sample 
from the moment it was inserted in the heated part of the instru¬ 
ment, the sintering percentage that was evaluated measuring the 
height variations of the sample, and the area variation that mea¬ 
sures the crossing area of the sample. The most important trend 
is that of the sintering percentile. A reduction under the 95% range 
is symptomatic of ashes sinterization and therefore slagging is 
imminent. The test outlined a temperature of 750 °C, that was low¬ 
er than the temperatures recorded during the tests. The presence 

Table 1 

Biomass proprieties. 


Proximate analysis (%wt) 

Volatile matter 

79.1 

Ultimate analysis (%wt) 
C 

47.7 


Fixed carbon 



Physical and chemical proprieties 
Chip equivalent diameter 


Biomass moisture 

Dry Biomass higher heating value 


0.0156 m 
0.4 kg/m 3 

0.63 
0.03 

15.67 MJ/kg 


of slags was then considered as additional proof that the reactor 
working on 15-s loading was characterized by the channeling phe¬ 
nomenon. In the channels, where the equivalence ratio is high, the 
temperature rose and the slags were created even if the thermo¬ 
couples were recording lower temperatures in other parts of the 
reactor. 


2.4. Mathematical modeling 


The mathematical model was adapted from two different liter¬ 
ature models with appropriate modifications and improvements. 
The first model was developed by Reed and Levie (1984); it pre¬ 
dicts the flaming pyrolysis zone length l p (m) and the char reduc¬ 
tion zone length l c (m) starting from biomass properties and 
gasifier dimensions with the following equations: 

Ip = Vftp (1) 

lc = V f t c (2) 


Vf = 


rh bi0 (A g F d (1 - F„)) 
A g 


(3) 


where V f is the fuel velocity (m s '), t p is the pyrolysis time (s), t c is 
the char reduction time (s), m bi0 is the input biomass flow (kg s '), 
A g is the area of the gasifier (m 2 ), F d is the density of the biomass 
(ad) and F v is the void fraction in the biomass (ad). The pyrolysis 
time t p was obtained by the following equation (Reed and Levie, 
1984): 


F d V(h p + F m h w ) 
Aq 


(4) 


where V is the volume of the biomass particle (m 3 ), A is the surface 
area of the biomass particle (m 2 ), h p is the heat per unit mass re¬ 
leased by the pyrolysis process (kj kg -1 ) at the temperature T s (K), 
h w is the water latent heat of vaporization (kj kg -1 ) and q is the heat 
transfer rate per unit area in the pyrolysis process by radiation 
(kW m 2 ). The values of h p and h w were tabulated by Reed and Levie 
(1984) starting from the data obtained by Huff (1982). The heat 
transfer rate q can be obtained with a weighted average calculation 
from the data reported by Reed and Levie (1984), where the weight 
of the calculation comes from the moisture of the biomass, the sur¬ 
face temperature and the surface area of the biomass particle. The 
char reduction time was calculated assuming a constant height H 
of the fixed bed (m) from the following equation: 


t c =(H-l p )/V f 


(5) 


Eqs. (l)-(5) allow a calculation of the bed and zone heights, the 
time of the pyrolysis process, the char reduction process assuming 
a constant pyrolysis surface temperature and a constant biomass 
consumption. 

An additional approach, presented by Wang et al. (1993), was 
exploited to estimate the syngas composition, the air inlet flow, 
the syngas output flow and the charcoal and water production. 
This model does not consider the tar production, the temperature 
of the reduction zone and that the oxidant equivalence ratio are 
fixed. It is based on the following general reaction of gasification 
(Wang et al., 1993): 


CH a O^ + y0 2 + zN 2 + wH 2 0 = <p\C+ </> 2 H 2 + 04 CO + ^> 4 H 2 0 

+ ^5 C0 2 + tj> 6 C H 4 + faN 2 (6) 


where the parameters a and /i can be calculated from the ultimate 
analysis of the biomass normalizing with the specific weight of car¬ 
bon, y and z are the moles of molecular hydrogen and nitrogen for 
mole of biomass, these values can be calculated from the reaction 
balance of the biomass oxidation considering the equivalence ratio 
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ER defined as the ratio between the air flow and the stoichiometric 
air flow (Reed and Das, 1988); w is the number of moles of water 
per mole of biomass evaluated from the biomass moisture F m , fa 
are the molar numbers of the products of the gasification per mole 
of biomass. Eq. (6) can be written in mass flow mode and it allows 
to calculate the syngas composition, the inlet and output flow rate 
assuming a constant biomass consumption. The variables fa are the 
unknowns of the problem. At the end of the flaming pyrolysis zone 
(subscript 0) these variables were calculated by solving the linear 
system below (Wang et al., 1993) assuming <j> 20 — 0, r/> 30 = 0 and 


07,0 = Z- 

01,0 + 05,0 + 06,0 = 1 

(7) 

204,o + 40 6,0 = a + 2 w 

(8) 

04,0 + 205,0 = 2 y + P + w 

(9) 

04,0 = 205,0 + W 

(10) 

where X = 1 is the ratio between water vapor and carbon dioxide 
formation (Wang et al., 1993). In the reduction zone the tempera¬ 
ture decreases from 1200 K to 1000 K. Here, char is consumed grad¬ 
ually due to char-gas reactions and gas-gas reactions as described 
by the following chemical reactions (Wang et al„ 1993): 

1. Boudouard reaction: C + C0 2 <=* 2CO 

2. Water-gas reaction: C + H 2 0 <=* H 2 + CO 

3. Methanation reaction: C + 2H 2 CH 4 

4. Steam reforming reaction: CH4 + H 2 0 <=> CO + 3H 2 


Surface reactions 1 and 2 involve single gas molecules, while two 
molecules are involved in surface reactions 3 and 4. The Lang- 
muir-Hinshelwood mechanism (Laidler, 1987) is applied to calcu¬ 
late the net rate for reactions 1 and 2 given by Eqs. (11) and (12), 
and the Langmuir-Rideal mechanism (Laidler, 1987) was adopted 
to calculate the net rate for reactions 3 and 4 given by Eqs. (13) 
and (14) as follows: 

,, , 05-0^/Wl) /01,oV /3 /0, A 

El 2 (ffi + l/p)0i V 0i 7 VM>y 

(11) 

r k 04 - 0203/CfW /01,O\ V3 / 0i \ 

2 “ 2 E,L 2 (K. + i/p)0i v 0. y ypdp) 

(12) 

t . , 0^-0 6 fVK P 3 (hA l/3 (4>i) 

3 03 P^ELCKf + VP)0i V 0i ) \Pdp) 

(13) 

„, ,, 0406 - 0302/(P^W fhflY /3 f 0i \ 

4 “ 4 EL(K. ' 1/P)0f \0fi VPV 

(14) 

-Wa 

#■ 

(15) 


where Vi is the net reaction (mol s 1 ), k ai is the apparent rate con¬ 
stant, K pi is the equilibrium constant and K, is the adsorption con¬ 
stant for the reaction i, p is the pressure in the gasifier (atm), P, ; , 
is the pressure related constant obtained as described in Eq. (15), 
fa 0 is the initial molar number of atomic carbon at the beginning 
of the reduction zone, p is the density of the atomic carbon 
(kg dm -3 ), d p is the initial equivalent diameter of the char particle 
(m). The equilibrium constants I< pi were calculated by the JANAF 
thermochemical tables (Sharma, 2008) and the adsorption con¬ 
stants Ki were taken from Everson et al. (2006) both at a fixed 


Table 2 

Model parameters. 

Parameter 

Flaming pyrolysis temperature T s 

Fixed bed height H 

Biomass coefficient a 

Biomass coefficient /i 

Equivalence ratio ER 

Reduction temperature T 

Pressure in the gasifier p 

Tar higher heating value HHVt, 


Value Unit 

1173.5 K 

0.6 m 

1.52 

0.67 

0.3 

1073.5 K 

1 atm 

24.62 MJ/kg 


Table 3 

Comparison between the model output and the 4-s-long loading experimental data. 


Variable 


Biomass consumption m M(1 519 

Syngas higher heating HHV S 4.99 

Syngas volumetric flow Vsyngas 112 


Air volumetric flow rate v mr 

Char flow rate m char 

Tar + water flow rate m tor+ „ 

Cold gas efficiency r/ gas 


5.00 

11.9 


Unit 


519 g/s 

5.25 MJ / 

Nm 3 

103 1/s 

91.2 1/s 

6.90 g/s 

12.7 g/s 

66.6 % 


reduction temperature T (K). The apparent rate constants were cal¬ 
culated by the following Arrhenius equation (Dasappa et al., 2011): 

k ai = A exp ( 16 ) 

where A, is the pre-exponential factor (s '), E ai is the activation en¬ 
ergy (kj mol-1) for reaction i, T is the reduction temperature and R 
is the universal gas constant (kj mol-1 K— 1). A, and E ai of each reac¬ 
tion were taken from the regression reported in Wang et al. (1993). 
The molar number of the products at the end of the reduction zone 
was obtained by integrating the following differential equations 
system from zero to t c : 


v x + v 2 + v 3 
-v 2 + 2v 3 -3v 3 
-2v\ - v 2 - v 4 

V2 + V4, 

V\ 

-v 3 + v A 


(17) 


Table 4 

Comparison between the model output and the 15-s-long loading experimental data. 


Parameter Value 

Channels area Ai 0.318 

ER of the channels area ER, 0.667 

Pyrolysis area A 2 0.350 

ER of the pyrolysis area ER 2 0.048 

Variable Model 

output 


Biomass consumption 
Syngas higher heating 

Syngas volumetric flow 

Air volumetric rate 
Char flow rate 
Tar + water flow rate 
Cold gas efficiency 



Experimental Unit 


531 g/s 

3.97 MJ/ 

Nm 3 

103 1/s 

91.2 1/s 

7.77 g/s 

21.8 g/s 

66.6 % 
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The molar flow rate of syngas, char, and water at the end of the 
reduction zone were evaluated using the molar mass of each prod¬ 
uct and assuming steady state conditions. Starting from the gas 
composition one can evaluate the syngas heating value HHVsy^as 
(MJ Nm~ 3 ) thanks to Waldheim and Nilsson equation (Waldheim 
and Nilsson, 2001). HHV Syngas is used to obtain the cold gas efficiency 
of the reactor: 


_ VsyngasHHVsyngas 

^ “ m bi0 HHV bi0 


(18) 


where V syngas is the syngas flow rate (Nm s 1 ). Moreover the higher 
heating value of the biomass HHV bi0 (MJ kg -1 ) was calculated by Eq. 
(19) starting from the moisture of the biomass and the higher heat¬ 
ing value of the dry biomass HHV hin d ^ (MJ kg -1 ) estimated experi¬ 
mentally by the Mahler bomb calorimeter (ASTM, 2007): 


HHV bi0 = HHV biodr y(\ - F m ) - 2.257 F m (19) 


The tar mass flow rate rh t „r (kg s 1 ) was evaluated using an energy 
balance equation: 


as ~ V sy ngas,expHHVsyngas,exp 

HHV [arexp 


(20) 


where V syngas , exp is the experimental volumetric flow rate, 
HHVsyngas,exp is the experimental syngas higher heating value and 
HHVtar,exp is the experimental tar heating value. With this simplified 


approach, the tar production is evaluated looking at the differences 
between the quality of the modeled gas and the real gas. The basic 
assumption behind this method assumes that the real gasifier is not 
able to completely convert the biomass chemical energy into the 
gas chemical energy. Therefore the only loss in the conversion pro¬ 
cess is the tar production, which represents a way to discharge en¬ 
ergy from the system. 

A Matlab™software was used to implement the mathematical 
model discussed below. The biomass adopted in this study was 
wood-chip of Popolus Alba (Tutin and Heywo, 1976) and Table 1 
collects its chemical, physical and geometrical proprieties. Table 2 
shows the model input parameters, instead Tables 3 and 4 summa¬ 
rize the comparison between experimental data and model output. 
The HHV hi0: dry and the HHV tar were evaluated by a calorimetric 
analysis with the Mahler bomb calorimeter (ASTM, 2007). 

A heuristic upgrading of the present model was employed to 
consider the presence of channels in the 15-s-long loading opera¬ 
tion mode. The channels were modeled considering them as a part 
of the reactor with a total area A d , while the rest of the reactor bed 
was assumed to have a section with total area A 2 — A - A u Biomass 
consumption in each sub-reactor was considered proportional to 
its area. Each part was simulated with different boundary condi¬ 
tions: ERj, mi - mA,/A for the first part; ER 2 , m 2 = mA 2 /A for 
the second part. A method that randomly changed the boundary 
conditions was applied to verify if the channeling theory could ex¬ 
plain both the high water and tar production. Table 4 collects the 


(a) 



(b) 



Fig. 3. Thermal stratigraphy in the 4-s-long loading test (a) and in the 15-s-long loading test (b). 
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boundary conditions that better fit with the experimental output. 
This approach went under the name of the ‘inception model’ be¬ 
cause it simulated the presence of a gasifier within a gasifier. 

3. Results and discussion 

3.1. Temperature 

The thermal stratigraphy of the top part of the gasifier is re¬ 
ported in Fig. 3. The temperatures were lower than what literature 
suggests in both the loading conditions. The temperature in the 
center of the reactor should have been higher than that which 
was attained by the thermocouples because the reactor diameter 
was big and the biomass could not properly reach the peripheral 
part of the gasifier. Here the biomass created a small ash annulus 
that protected the thermocouples from the higher temperature of 
the gasifier. During the plant maintenance operations, the ash layer 
was measured and it proved to be a maximum of 3 cm thick. The 
reduced thickness justifies that this phenomenon was in fact ig¬ 
nored in the model. 

Fig. 3 shows the different temperature trends at the top of the 
reactor for the two loading conditions. In the 15-s-load feed 
(Fig. 3 (b)), the large amount of biomass for every “batch” de¬ 
creased the reactor temperature due to the sensitive heat absorbed 
by the biomass and choked the flaming combustion reactions on 
the top surface of the reactor bed. The 15-s-long loading test 
showed lower temperatures and small fluctuations, especially for 


the thermocouple placed under the flaming pyrolysis zone (TC- 
3). In fact, the reactor was not able to recover the high temperature 
suggested by literature for these zones. The 4-s-long loading con¬ 
dition created greatly fluctuating high temperatures, even for the 
deeper thermocouple. These variations were related to an irregular 
gasification zone: the width of the peaks was about 60 s, this was 
the stand-by period for every 5 s run of the stirrer. In the 15-s test, 
the thermal substratum had a low temperature and was quite 
homogeneous: for this reason the stirrer did not create significant 
variations, but during the 4-s test (Fig. 3(a)) this substratum was 
characterized by higher temperatures ranging from 500 to 
700 °C, which are the extremes of gasification reaction tempera¬ 
tures as reported by Reed and Das (1988), Higman and van der 
Burgt (2008) and Knoef (2005). 

Such a non-homogeneous zone was greatly modified by the 
passage of the stirrer, thus justifying the significant variations. 
Fig. 4 shows two trends obtained acquiring the superficial temper¬ 
ature of the reacting bed with the IR thermometer. The thermom¬ 
eter was a dual laser infrared thermometer with focus 50:1 inches 
fixed on the air inlet chimney, focusing it on the top of the reacting 
surface. 

The differences between Fig. 4(a) and (b) consist in the maxi¬ 
mum temperature reached by the surfaces; furthermore the tem¬ 
perature recovery was slower in the 15 s loading. The 15 s 
loading showed 9 big peaks during the 1000 s of the screenshot, 
it resulted in a 15 s long loading every Ills. For the 4-s test this 
value decreased to 77 s on average. In order to estimate the 


(a) 



(b) 



Fig. 4. Superficial temperature in the 4-s-long 
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T [°C] Sintering [%] Poplar wood, chips, klin closed at 5,00 °C 



amount of time in which the combustion reactions were locally 
choked due to the sudden loading, the temperature trends were di¬ 
vided into two zones: above and below 250 °C. This value was the 
average autoignition temperature for wood (Babrauskas, 2001 ); for 
this reason, during the periods with a reactor surface temperature 
above 250 °C, a wood parcel on the top surface could start its flam¬ 
ing pyrolysis if it was able to reach the right amount of oxidant. Re¬ 
sults of this analysis are reported in Table 5. The 4-s loading 
condition resulted in a surface reactor temperature above the 
autoignition one for the 67% of the duration of the test. During 
the 15-s-loading campaign this value was lower (39% of the dura¬ 
tion of the test). A lower value can be related to a higher tar pro¬ 
duction because temperatures under 250 °C, associated with the 
choking effect of long loading, generate a more similar reaction 
to the torrefaction than to that of flaming pyrolysis (Prins et al., 
2006). 


3.2. Power output and propeller frequencies 

The power plant studied here is controlled on a power output 
basis. Once the power output was set, the controller continuously 
changed the impeller frequencies, increasing or decreasing the 
amount of syngas drawn from the reactor and sent to the engines, 
in order to maintain the fixed power generation. Changing the syn¬ 
gas flow, constantly modified the reaction conditions. This phe¬ 
nomenon created oscillations in the propeller frequency and in 
the instantaneous power output registered from the controller 
data logger. 


Table 5 

Superficial temperature of the reacting bed comparison. 

Variable 4 s load 15 s load Unit 

Fraction time over 250 °C t hM 67.89 38.6 % 

Minimum average temperature T mi „ 42.6 31.8 °C 

Maximum average temperature T mm 896 739 °C 

Average temperature T mM „ 422 259 °C 


In order to quantify differences in the instability of the whole 
system under the two loading conditions, the instantaneous values 
of the impeller spinning frequencies and power output were ac¬ 
quired. The mean value and the standard deviation were calculated 
and reported in Table 6. For the 15-s loading the impellers spun 
faster on average. This can be related to a poorer quality of syngas 
because a greater amount of gas was necessary for the same power 
output. Moreover, the standard deviation of the spinning fre¬ 
quency was also higher in this case. This suggests that the great 
thermal and physical oscillations at the top of the reactor influ¬ 
enced the behavior of the propellers too. The result obtained for 
the power output was different and Table 6 shows a lower stan¬ 
dard deviation of this parameter in the 15-s loading. This phenom¬ 
enon can be explained looking at the temperature trends in Fig. 4: 
during this loading condition the temperatures rose slowly after 
every reactor refill. For this reason, even if the impellers were 
forced to greatly change their spinning velocities due to the large 
changes in the reacting conditions, these changes occured slowly, 
so the impellers could properly “follow” the changes maintaining 
a more constant power output. 
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Propeller frequencies comparison and engines output comparison. 

Variable 4 s load 15 s load Unit 

Minimum average f mln 

propeller 
frequencies 

Maximum average f max 

propeller 
frequencies 

Average propeller f mi , an 

frequencies 

Standard deviation of std,/ meon 

the propeller 
frequencies 

Minimum average P min 

electrical power 

Maximum average P max 

electrical power 

Average electrical power P mea „ 

Standard deviation of std, P mem 

the electrical power 


3.3. Model and experimental output 

Tables 3 and 4 collect the main outputs for the two different 
loading conditions compared with the analogous simulated sys¬ 
tem. The more important data for channeling monitoring were 
the water production and the tar production, they can be con¬ 
sidered in contrast with each other because a higher equiva¬ 
lence ratio causes low tar and high water production due to 
the preponderance of the combustion reactions. The opposite 
condition is related to a low equivalence ratio that produces oils 
and tars instead of water. The 15-s-loading campaign resulted in 
a high production of both water and tar. The scientific explana¬ 
tion came from the ‘inception’ model of the reactor, which sim¬ 
ulated a reactor characterized by the coexistence of combustion 
paths in channels while the rest of the reactor works with a low 
equivalence ratio condition. Table 3 shows the effectiveness of 
the model adopted for the prediction of the main parameters 
that characterized the reactor without channeling, while Table 4 
shows how the model divided the reactor in order to match the 
experimental data. The model suggested an almost equal divi¬ 
sion of the reactor area with completely different equivalence 
ratios. Under these conditions half of the gasifier acted as a 
pyrolysis reactor producing high amounts of tar. The rest of 
the gasifier was responsible for the water and ash slags 
production. 

Despite the presence of a double bladed stirrer in the reactor, 
which had the purpose of increasing the fuel flow and reducing 
the channeling (Patil et al., 2011; Ma et al., 2012), this phenom¬ 
enon still affected the gasifier. This could have been caused by 
a peripheral location of channels, but the low temperatures re¬ 
corded through the reactor jacket suggested abandoning this 
hypothesis. A possible interpretation of the presence of the chan¬ 
neling even in a gasifier equipped with a stirrer consisted in a 
synergy of two factors: first of all, due to the particular air inlet 
design, the flaming pyrolysis zone was not flat and the stirrer 
could contribute to damaging the fragile equilibrium of this ther¬ 
mal zone. Moreover the channeling formation rate may have been 
faster than the capability of the stirrer to desegregate the chan¬ 
nels. These observations outlined the importance of working on 
control parameters and modeling before designing structural 
modifications to the plant or the reactor; the results reported 
here show how the subdivision of the loading in small amounts 
drastically reduced the tar formation and the water production 
in the reactor. 


3.4. Biomass consumption and char production 

Comparing the two tables, the amount of biomass consumed 
per hour did not seem to be influenced by the loading conditions 
and it depended only on the set power output of the whole system. 

The charcoal production required particular careful because it 
was similar under the two conditions. This is somehow an unex¬ 
pected result, because a higher amount of energy was presumed 
to be discharged by the char in the 15-s-loading case. For this rea¬ 
son, the different charcoals were tested in a bomb calorimeter in 
order to estimate their high heating value. The charcoal produced 
during the 15-s loading test had a HHV of 24.87 MJ kg -1 , higher 
than the 17.89 MJ kg 1 obtained testing the charcoal produced 
during the 4-s loading test. This phenomenon was caused by the 
char discharging system, which worked maintaining a fixed pres¬ 
sure drop at the bottom of the reactor, ignoring the fact that the 
char had or had not finished reacting. 

4. Conclusions 

The experimental campaign showed a significant relationship 
between loading parameters and channeling phenomenon in 
downdraft gasifiers. Contemporaneously to these investigations, 
the channeling produced by loading the reactor with 15 s-long 
batches found a theoretical counterpart thanks to the “inception” 
mathematical model. This model explained the effect of large loads 
that choke the superficial reactions, forcing the air to run through a 
few paths in the reacting bed. Channels generated in such a way 
are characterized by a high equivalence ratio where most of the 
water is produced, while the rest of the reactor runs with a low 
equivalence ratio, producing a great amount of tar. 
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21.4 21.3 Hz 

37.8 38.9 Hz 

29.9 30.78 Hz 

2.71 3.08 Hz 

119.5 116.6 kW 

207.7 204.6 kW 
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4.79 4.17 kW 








